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A number of spectacular experimental 
anomalies 1 ' 2 have recently been discovered 
in certain cuprates, notably L&2 - x Ba x Cu04 and 
Lai.6- a; Ndo.4Sr a; Cu04, which exhibit unidirec- 
tional spin and charge order (known as "stripe 
order"). We have recently proposed to inter- 
pret these observations as evidence for a novel 
"striped superconducting" state, in which the 
superconducting order parameter is modulated 
in space, such that its average is precisely zero. 
Here, we show that thermal melting of the 
striped superconducting state can lead to a 
number of unusual phases, of which the most 
novel is a charge 4e superconducting state, with 
a corresponding fractional flux quantum hc/Ae. 
These are never-before observed states of matter, 
and ones, moreover, that cannot arise from the 
conventional Bardeen-Cooper-Schrieffer (BCS) 
mechanism. Thus, direct confirmation of their 
existence, even in a small subset of the cuprates, 
could have much broader implications for our 
understanding of high temperature supercon- 
ductivity. We propose experiments to observe 
fractional flux quantization, which thereby could 
confirm the existence of these states. 

It is widely accepted that the mechanism of supercon- 
ductivity in the cuprate high temperature superconduc- 
tors is different than in conventional superconductors, 
and that the "normal state" above T c is anything but 
normal. In contrast, the low energy low temperature 
properties of the superconducting state appears to be 
largely accounted for by the conventional BCS theory, 
albeit with a d-wave order parameter and a much sup- 
pressed superfluid density. In short, the perception is 
that the superconducting (SC) phase is simple while the 
"normal" phase is not. However, recent experiments^ 
in the "stripe ordered" materials, in which a novel form 
of highly two-dimensional partial SC order was found, 
suggest that the SC phases may have unconventional as- 
pects. We have interpreted^^ these anomalies as evi- 
dence of a new type of order, "striped superconductivity" 
or unidirectional pair density wave (PDW) order. Sub- 
sequent experiments^ have provided additional indirect 
evidence in favor of this interpretation. 

In the present paper, we study the thermal melting 
of a striped superconducting groundstate by the prolif- 
eration of topological defects leads to a complex phase 
diagram, shown in FigfTJ with several interesting phases, 
of which the most novel is a charge 4e superconducting 



state. The existence of such a state can be directly estab- 
lished experimentally by observing magnetic flux quan- 
tization with period hc/Ae, half that of the usual super- 
conducting flux quantum, in a variety of experimental ge- 
ometries, as shown in Figf2] Moreover, it is likely that the 
charge 4e SC, resulting from partial melting the striped 
SC, has a finite density of states for gapless quasiparti- 
cles, reflecting the fact that it does not arise from the 
Bose condensation of literal four electron bound states. 

In related developments, Agterberg and Tsunetsugu&2 
and Radzihovsky and Vishwanathfi^ have discussed re- 
lated phenomena in connection with Larkin-Ovchimiikov 
(LO) phases in, respectively, the heavy fermion mate- 
rial CeCoIns at high fields and in partially spin-polarized 
ultra-cold atomic gases. The states we discuss here dif- 
fer from the LO states in that there is no explicit time- 
reversal symmetry breaking. Also the structure of the 
striped SC reflects the discrete point group symmetry of 
the underlying (cuprate) lattice. 

A striped superconductor is a unidirectional PDW 
state in which a spin singlet superconducting order pa- 
rameter, A(r) oscillates in space with a wave vector Q, 11 
(shown schematically in Fig. [I}: 

A(r) = A Q (r) e iQr + A_ Q (r) e^ * . (1) 

Of course, there will always be higher harmonics of the 
order parameter at wave vectors nQ, but in the ordered 
state, these are slaved to the fundamentals and hence are 
not independent dynamical degrees of freedom. Gener- 
ally, therefore, we will not treat them explicitly. How- 
ever, two subsidiary order parameters play a special role 
in the thermal melting: Charge density wave (CDW) or- 
der, p(r) — po + /9k (r) e lK r + p K (r) e _lK r , appears as a 
second harmonic of the fundamental ordering, with wave 
vector K = 2Q, where /3k (r) °c A*g(r)Aq(r). Sim- 
ilarly, charge 4e SC order, represented by the complex 
scalar field A4 e (r) oc A_Q(r)Acj(r), occurs parasitically 
along with the fundamental striped SC order. 

Although the stripe ordered superconductor need not 
be associated with a non-zero magnetization, it is in 
many ways similar to a LO state.— As such it has a 
U(l) x U(l) symmetry corresponding to the independent 
uniform shifts of the phases of the complex components 
of the order parameter: 

A±q (r) -, e '^A ±Q (r)=e'( tt «A ±Q (r), (2) 
PK(r) ^e l20 /3 K (r) with = (0 Q - 0_ Q )/2 
A 4e (r) ^e l2e A 4e (r) with d=(d Q + 9_ Q )/2 
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Much as other stripe electronic liquid crystal states, the 
PDW also typically breaks the point group symmetry; 
to be concrete, we will consider the case of a tetragonal 
crystal, in which the choice of direction for the unidirec- 
tional order breaks the C 4 . u point group symmetry down 
to C-2. In analogy with classical liquid crystals, wc refer 
to this as nematic ordering, although since the symme- 
try breaking is in fact, discrete, it is an "Ising nematic." 
This is the minimal set of broken symmetries associated 
with the striped superconductor. For simplicity, we will 
consider only the interplay of charge stripe and SC stripe 
orders, and ignore the striped spin order, even though all 
three are intertwined where striped superconductivity is 
conjectured to occur in the cuprates^^ 

We consider a system which has a striped supercon- 
ducting groundstate, and consider the manner in which 
thermal fluctuations gradually restore all the broken sym- 
metries as the temperature, T, increases. Deep in the 
PDW phase, other than in a vortex core, we can ig- 
nore fluctuations in the magnitude of the order param- 
eter, and write A ± q(r) = A sc exp [i (8(r) ± 4>(r)) }, 
Pk(t) = Pifexp [i2(/>(r)], and A 4e (r) = A 4e exp \i26[r)]. 
Because the phase of the CDW order is subject to pinning 
by quenched disorder, the PDW is considerably more 
fragile than a uniform SC state, and quenched random- 
ness easily leads to an XY superconducting glass phase. 
(Unless stated otherwise, below we ignore the effects of 
disorder.) Moreover, the striped superconductor state 
admits novel topological excitations, which play a cen- 
tral role in the thermal melting of the groundstate. 

Deep in the striped SC phase, treating the system as 
2D, and taking Q in the x direction, the effective Hamil- 
tonian for the low energy thermal fluctuations is 

K[0M = yWl' + f (Dovf, (3) 

D s = a s [~ihd x - (2e/c)A x ]x + a^l-ihdy - (2e/c)A y ]y, 
D c = a c [— ihd x ]x + a~ 1 [—ihdy]y, p s and k are, re- 
spectively, the superfluid stiffness and the CDW elastic 
constant, and a s and a c are the corresponding (finite) 
anisotropics.— We have shown, explicitly, the coupling 
to an external gauge field, but henceforth, unless other- 
wise specified, we will take A = 14 . The correlation 
functions of the PDW, CDW and charge 4e SC order 
parameters in this phase are 

(A ±Q (r)A ±Q (rT) cx K s (r - r')K c (r - r') 

<P2 Q (r)p 2Q (r')*)cx[X c (r-r')] 4 

<A 4e (r)A 4e (r')*}oc[^ s (r-r')] 4 (4) 

where 

K s (r)~ [(xa s ) 2 + (a-'y) 2 ]-^ 

K c {r)~[{xa c f + {a- l yf}-^ 

r/ s = 2-wT/ps r/ c = 2-kT / n. (5) 

Since the order parameters A±q(i") must be separately 
single valued, their phase fields must be invariant un- 
der the transformations 0±q(t) — > 9±q(t) + 27rm±Q, 



where m±Q are integers. Correspondingly, the fields 
9(r) and ip(r) must obey the conditions 6(r) — > 9(r) + 
7t(toq + m_ Q ), iys(r) — > (f(v) + 7t(toq - m_ Q ). The 
integers to±q then classify the topological excitations 
supported by the PDW state: vortices with topological 
charge q s = (toq + m_Q)/2 and dislocations with topo- 
logical charge q c — (jtiq — to_q)/2. We have three types 
of topological excitations (q s , q c ): 

1. Full vortices with q s = ±1 and q c = 0. 

2. Double dislocations, with q c = ±1 and q s = 0. 

3. Half- vortices, q s = ±1/2, bound to single disloca- 
tions, q c = ±1/2. 

Much as in the case of the well understood theory 
of the Kosterlitz-Thouless (KT) phase transition ) 15 ! 16 ^ 
the thermodynamic behavior of this system is also rep- 
resented by a generalized (neutral) Coulomb gas of vor- 
tices {q s } and dislocations {<z c }, both with logarithmic 
interactions, where the strength of the interactions be- 
tween vortices (and anti- vortices) , both fractional and in- 
tegral, is controlled by the superfluid density p s and the 
interaction between dislocations (and anti-dislocations) 
is controlled by the CDW stiffness k. Let us denote by 
9(is,qa) t ne fugacity of a topological excitation with topo- 
logical charges {q s ,q c ). In the dilute limit, in which the 
fugacities are small, the partition function of the gen- 
eralized vector Coulomb gas 1 ^ can be represented by a 
sine-Gordon type effective field theory (for a review see 
Refs»i£), that in this case requires two fields 9 and ip, 
the dual of the superconducting phase 9 and the CDW 
phase (p. The effective Hamiltonian density of the 2D 
field theory, dual to the degrees of freedom of Eq. §5§ , is 

2p s \ / 2k 
-5(i,o) cos(2tt(?) - g (0A) cos(2tt^) 
-2.9 ( i , i) cos(7rf?) cos(tt^) (6) 

There are three distinct pathways for the thermal melt- 
ing of the PDW stated by condensation of its topolog- 
ical excitations, as shown in the schematic phase dia- 
gram in FigfT] In constructing this figure, the location of 
the phase boundaries were estimated by computing the 
scaling dimensions of the cosine operators at the non- 
interacting fixed point, A gsj ^ c = i^(q 2 p s + q 2 n), and then 
identifying the lines, A gs;9c = 2, at which each of the 
cosine operators in Eq.© first becomes relevant. Alter- 
natively, a mean-field phase diagram can be obtained by 
treating Eq.© in a self-consistent phonon approxima- 
tion. More sophisticated RG treatments of a Coulomb 
gas system with precisely the same formal structure as 
the above have been discussed previously for a model 
of the thermal melting of a charge-spin stripe state^S 
and the context of spinor S = 1 Bose condensates^ 
All approaches yield the same topology of the phase dia- 
gram as the one shown in Fig[TJ although there are some 
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FIG. 1: Schematic phase diagram: KT transitions between 
the PDW, CDW (stripe), nematic, and charge 4e SC, and the 
(Ising) nematic-isotropic transition. Mi and M2 are multi- 
critical points (see text). PDW, CDW, and charge 4e SC 
orders have QLRO in their respective phases, and are short 
ranged in the nematic phase. 



unresolved differences concerning the shape of the vari- 
ous phase boundaries and the nature of the multicritical 
points, Mi and M 2 . Moreover, since the specific prob- 
lems being addressed in^l and^S were different than those 
being considered here, the interpretation and the physics 
of the phases was different. 

a. PDW Phase: At low temperatures, all three co- 
sine operators are irrelevant and the topological excita- 
tions are uncondensed. The system is in a PDW state 
with quasi-long range order (QLRO), as given in Eq. (J4|. 

b. Charge 4e SC Phase: This phase is accessed when 
the double dislocations proliferate, i.e. when the opera- 
tor cos(27T(^) becomes relevant, Hence, the PDW-Charge 
4eSC phase boundary is the line A ,i = 4nn/T = 2, 
or — = 27t^ (see FigQ]) In the charge 4e SC phase, 
the field (p is pinned at integer values n v G Z, and has 
massive (gapped) fluctuations. In this phase, the dislo- 
cations of the CDW are screened, and the dual the CDW 
phase field Lp has wild fluctuations, leading to exponen- 
tially decaying correlations of both the PDW and CDW 
order parameters. In contrast, the field 9 remains gap- 
less, and the correlation function of the charge 4e SC have 
power-law correlations, with the same exponent rj4 e . The 
phase transition from the charge 4e SC to the disordered 
non-superconducting nematic higher temperature phase 
proceeds through the subsequent unbinding and prolifer- 
ation of vortices with fractional (half) topological charge. 

c. Stripe (CDW) Phase: This phase is accessed 
through the proliferation of SC vortices with integer 
topological charge when the operator cos(27r6' + ) becomes 
marginal, A 1>0 = np s /T = 2. There is CDW QLRO, 
with correlation functions that decay with a power law 



with an exponent r\CDW and no SC order of any type. 
In this phase the field 9 is pinned at the integer values 
n$ G Z, and its fluctuations are massive. The phase 
transition from this phase to the normal (nematic) phase 
occurs by the unbinding of single dislocations whose frac- 
tional vortex charge is screened in this phase. 

d. Direct PDW -nematic normal state phase transi- 
tion: This transition proceeds through the liberation 
of half-vortex-integer dislocation composite excitations, 
i.e. when the topological excitations (±5, ±5) condense. 
This phase boundary is the M±M 2 line shown in Fig. Q] 
determined by the condition A ± i ± 1 = 2. 

e. Multicritical Points: The schematic phase dia- 
gram of FigQ] has two multi-critical points, M\ and 
M2, whose existence follows from the topology of the 
phase diagram. From the structure of Eq. [5J it also 
follows that there exist two special critical points with 
an emergent higher symmetry where two operators be- 
come marginal simultaneously: (0, 1) and (|,±|) at M[, 
and (1,0) and (^,±i) at M 2 . At these points, the gen- 
eral U(l) x U(l) symmetry of the model is enlarged to 
an SU(3)i symmetry— and in their vicinity the cor- 
relation length exhibits KT behavior with a modified 
exponent . 1 20 The simple scaling arguments used to con- 
struct FigQ] and the RG analysis ir*2£ suggest that the 
high symmetry and multicritical points are one and the 
same, Mj — Mj, while the self-consistent phonon cal- 
culation and the RG treatment of^i suggest that they 
are distinct, with M[ and M' 2 lying part way along the 
segment of the phase boundary between M\ and M 2 . 

While the present discussion of thermal melting was 
confined to 2D, it can be readily extended to the case of 
a 3D layered material with sufficiently weak inter-layer 
coupling. Some aspects of this extension depend on de- 
tails of the interlayer geometry^. For instance, depend- 
ing on whether or not the inter-layer coupling frustrates 
the stripe ordering, the low temperature state may or 
may not have a gentle, time-reversal symmetry breaking 
superconducting spiral superimposed on the basic single- 
plane stripe order. But, in general, all of the phases in- 
dicated in FigQ] persist in 3D, with the QLRO replaced 
by true long-range order, and a crossover, very near the 
phase boundaries, to 3D criticality. 

At T = the power-law correlations of the PDW 
phase are replaced by true long range order and gap- 
less collective (Goldstone) modes. A simple mean-field 
treatment^ 2 ^ of the spectrum of Bogoliubov quasiparti- 
cles in this phase reveals that, in common with a CDW, 
it generically has an only partially gapped Fermi sur- 
face, and hence a finite density of zero energy states. In 
principle, it should also be possible for the PDW phase 
to undergo quantum melting, possibly also by a quan- 
tum analog of the vortex unbinding mechanism described 
abovei 2 ^ This is a more complex problem than the ther- 
mal melting we have discussed here. In particular, the 
quantum phase transitions are affected by the existence of 
gapless quasiparticles. Hence, the structure and quantum 
dynamics of the vortices (both integer and fractional) is 
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(a) (b) 

FIG. 2: Two schematic phase-sensitive experiments: a) a 
SQUID loop with a charge 4e SC, and b) a charge 2e striped 
SC-charge 2e uniform SC Josephson junction (see text). 



expected to be generally damped and anisotropic, which 
affects the physics of quantum melting. 

There are many unusual physical consequences of the 
nature of the PDW phase and of its thermally melted 
daughter phases. We have previously discussed^^ some 
remarkable bulk features of a fully ordered PDW phase, 
including the possibility of dynamical layer decoupling 
(in a 3D layered material), an anomalous sensitivity to 
quenched disorder, the existence (in principle testable by 
STM and scanning SQUID) of pinned half- vortices as- 
sociated with each disorder induced dislocation, and a 
tendency to formation of a superconducting glass phase 
which likely spontaneously breaks time reversal symme- 
try. Here we focus on "phase sensitive" measurements, 
which have not yet been explored, and offer the possi- 
bility of directly establishing the existence of either the 
fully ordered phase or its uniform, translationally invari- 
ant charge 4e SC descendant. 

Clearly, a superconducting SQUID loop, made of a 
charge 4e SC, will exhibit all the familiar features of a 
SQUID, but with a distinct half- flux quantum, hc/Ae re- 
placing the usual supeconducting flux quantum, hc/2e. 
From a practical point of view, however, it is undoubtedly 
easier to fabricate a SQUID loop in which only a single 
link consists of a charge 4e SC, as shown schematically 
in Figj2^. Such a SQUID will exhibit the same half-flux 
quantum as if it were an entire loop of charge 4e SC. In 
addition to the practical advantage, such a geometry will 
detect PDW related phases under a wider range of cir- 
cumstances. For instance, in the disordered phase near 
criticality, so long as the width of the putative charge 4e 
link is not large compared to the SC coherence length, it 
will simply act as a Josephson weak link in a charge 4e 
SC SQUID. Moreover, under many circumstances, even 
if the link consists of a PDW state, the SQUID will still 
exhibit hc/Ae periodicity. 

To see this, consider a Josephson junction between a 
striped SC and an ordinary SC, as shown schematically in 
FiglSb- The dependence of the free energy on the phase 



difference, A9, between the two SC can be expanded as 

oo 

F(A6) = S £ J n cos[nA6>] , (7) 

where S is the junction area, and J n involves tunneling 
processes of 2n electrons, so that generally J n decays ex- 
ponentially with increasing n. For the geometry shown 
in Fig. 2b, which is representative of the putative striped 
SC phase in L^-zBa^CuC^ near x — 1/8, J\ (and in- 
deed, all J2n+i) vanish identically, since a it phase change 
of the striped SC is equivalent to a translation by half a 
period. As a result, the (resistively shunted) Josephson 
coupling between the two SC is dominated by the intrin- 
sically smaller higher order coupling, Ji*^ Not only does 
that mean that a striped SC in the place of the charge 
4e SC link in FigJ^Ji, will act in the same way, it also 
means that all the standard characteristics of a single 
Josephson junction, including Josephson oscillations and 
Shapiro steps, will occur with twice the usual frequencies. 

J2 is proportional to J 2 , where J is the "bare" (micro- 
scopic) Josephson coupling. A quantitative estimate of 
J2 is difficult, since it depends exponentially on micro- 
scopic parameters. However, assuming that the coher- 
ence length £ is smaller than the period of the striped 
superconductor, an estimate based on an effective x — y 
model (similar to that of Reb£, but adapted for a c-axis 
junction between a striped and a uniform superconduc- 
tor) gives that J 2 ~ J 2 / | J'|, where J' is the inter-stripe 
Josephson coupling in the striped superconductor. Thus 
the energy denominator is the smallest energy scale of 
the single-plane problem. Moreover, since the total cou- 
pling energy is proportional to the area of the junction, 
it should be possible to use a large enough junction such 
that SJ2 is substantial (i.e., much larger than ksT). 

It should be possible to perform the proposed ex- 
periments using small crystals of La2- 2; Ba 2 ;Cu04 with 
x = 1/8, where there is already strong circumstantial 
evidence of the existence of a striped SC phase. The ob- 
servation either of hc/Ae flux quantization in a SQUID of 
the sort shown in Fig. 2a, or a charge 4e in the Josephson 
relation in a Josephson junction of the sort shown in Fig. 
2b, would constitute dramatic and direct evidence of the 
existence of these exotic superconducting phases. 

E. Berg, berez@stanford.edu, is the corresponding au- 
thor. 
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